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Abstract. We present measurements of the Goos-Hanchen shift in a two-dimensional 
dielectric microwave cavity. Microwave beams are generated by a suitable superposition 
of the spherical waves generated by an array of antennas; the resulting beams are then 
reflected at a planar interface. By measuring the electric field including its phase, 
Poynting vectors of the incoming and reflected beams can be extracted, which in turn 
are used to find the incoming angle and the positions where the beam hits the interface 
and where it is reflected. These positions directly yield the Goos-Hanchen shift. The 
results are compared to the classical Artmann result and a numerical calculation using 
Gaussian beams. 
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1. Introduction 



The Goos-Hanchen shift (GHS), the lateral shift of a totally reflected light beam at a 
dielectric interface, was discovered by Goos and Hanchen in 1947 pQ. It has subsequently 
been observed not only in optics, but also in acoustics [TJ|2] and neutron beams [3]. The 
quantum analog has been observed at graphene interfaces [I] . A simple theoretical model 
is due to Artmann [5 J ; there, the shift is calculated analytically for a beam consisting of 
two plane waves with slightly different propagation directions. Analytical calculations 
of the GHS for Gaussian beams have also been performed [21 [6] as well as numerical 
ones [TJ. 
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The GHS is difficult to measure directly for optical wavelengths because the shift 
is on the order of the wavelength and thus very small. Goos and Hanchen used multiple 
reflections; recently, a measurement technique [8] has been proposed which only needs 
one reflection and measures the GHS at a dielectric interface relative to a metallic one 
with a CCD camera. 

Because the GHS scales with the wavelength A of the beam, and thus vanishes 
in the limit A — > 0, it can be seen as a first-order wave correction to geometrical 
(ray) optics. Such corrections have recently been incorporated into a ray model for 
optical microresonators [9j [TOj El [T2J [13] in order to explain deviations of calculations of 
optical modes from ray predictions. The resulting "extended ray models" are interesting 
tools which can help to understand ray-wave correspondence in optical microcavities. 
Moreover, as microcavities can be viewed as open quantum billiards [Hj , this is directly 
related to the study of quantum-classical correspondence in open quantum systems [15] . 

Unfortunately, there are no experiments which directly measure the GHS in optical 
microcavities, as the measurement of modes with high spatial resolution is difficult in 
these systems. Moreover, all calculations have assumed reflection at planar interfaces 
or circular interfaces, even though the interfaces in microcavities typically have non- 
constant curvature. 

High-resolution measurements of electric fields, including the phase of the field, 
can be done easily in microwave resonators [HU [TO [18] ; they thus could become model 
systems where predictions for the influence of the GHS on modes in optical microcavities 
could be tested. It is also possible to do microwave GHS measurements at curved 
interfaces and directly study deviations from the results for planar or circular interfaces. 
The influence of boundary curvature is of particular interest in optical microcavities, 
as cavities with dimensions on the order of the wavelength of the optical modes have 
recently been fabricated [19]; in such systems, the boundary curvature is also on the 
order of the wavelength and typically not negligible. 

In this paper, we demonstrate that such measurements are in principle possible 
by considering the simplest scenario. We prepare microwave beams with a well-defined 
propagation direction and spatial extension, and measure their GHS upon reflection at 
a planar interface, comparing the results to both the well-known Artmann result and a 
Gaussian beam calculation. 

The paper is structured as follows. The experimental setup is described in section [2j 
and the generation of suitable microwave beams is described in section [31 The resulting 
beams and the GHS extracted from them are shown in section H] and compared to 
theoretical predictions. 

2. Experimental setup 

Experiments are done using a rectangular (50 cm x 100 cm) cavity made of Teflon with 
a refractive index n = 1.44. The vertical dimension of the cavity is h — 0.5 cm; modes 
with frequencies below z/ max = c/(2nh) 20.8 GHz with the vacuum speed of light c 
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can thus be regarded as two-dimensional. The electric field is measured using an Agilent 
8720ES vector network analyzer. The generated beam travels to one boundary of the 
cavity and are reflected; the GHS As can directly be extracted as the shift between 
incoming and outgoing (reflected) beams. Figure [I] shows the idea of GHS extraction 
from measured data. It is possible to extract Poynting vectors Si n and S out of the 
incoming and outgoing beams at each point r = (x,y) in space from the measured 
electric fields. We define the propagation direction of the incoming and outgoing beams 
as the average Poynting vectors (Si n ), (S out ). Because our measurements are done on a 
quadratic grid, they are given by 

N X Ny 



(Si 



(out)) = tt^tt 5^5Z Sin (° ut )^'^')' W 



y 1 = 1 J = l 



where the Xj, yj are the positions on which S in ( out ) values have been measured. N x ( y ) is 
the number of x(y) positions measured. 

The average Poynting vectors define the propagation direction of the incoming and 
outgoing beams. Together with the position r in = (x in ,y in ) on the incoming beam (r in 
is the average over all positions at which the incoming beam is measured), (Si n ) defines 
a straight line 

':H::) +( < s '"> < 2 » 

parametrised by t. The intersection of this straight line with the bottom of the Teflon 
plate y = yo yields the position x\ where the incoming beam hits the cavity boundary. 
Analogously, the position x<i where the outgoing beam starts at the boundary can be 
calculated by intersecting the straight line 

i)=fc) +t <*-> (3 » 

with y = yo- The GHS is then given by As = X2 — x\. 
3. Beam generation 

Microwave antennas produce spherical waves; for the GHS measurements, beams 
approximating plane waves have to be generated from these spherical waves. The 
generation of plane waves can be easily done because the superposition of N spherical 
waves ipj with wave number k and centers (xj,yj) on a straight line creates a beam 

3=N /2 j=N /2 exp Uk^/( X - xj) 2 + (y- yjf 

„> = E 2 ^ y) - E 2 ^ X - Xj r + (S - Vl r (4) 

which approaches a plane wave if N goes to infinity. A similar generation of a suitable 
beam from multiple spherical waves generated by microwave antennas has been used in 

EDI. 
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Figure 1. Extraction of the GHS As from the measured incoming and outgoing 
beams. (S; n ) and (S ou t) are averaged Poynting vectors of the incoming and outgoing 
beams and (x- m ,yi n ) and (x ou t,y ou t) position averages. x\ and X2 are intersections of 
the straight lines defined by the Poynting vectors and averaged positions with y = y$. 



The propagation direction of a plane wave generated from spherical waves can be 
influenced by adding a phase factor <j>(J) to each spherical wave. In this case, the 
generated beam is given by 

j=N/2 

il>(x,y)= ^yy m (5) 

j=-N/2 

By varying <f)(j), one can achieve propagation directions which lead to reflection under 
different angles of incidence. I.e., 4>(j) = j leads to reflection with an angle below the 
critical angle Xcr ~ 44 degrees for Teflon with n = 1.44, <p(j) = —j leads to reflection 
with an angle above the critical angle, and for <fr(j) = the reflection happens close 
to the critical angle. Unfortunately, there is no analytical formula for the relation of 
4>(j) and the resulting propagation direction; the choice of the <p(j) functions therefore 
remains somewhat arbitrary. 

A plane wave, however, does not experience the GHS upon reflection, as the GH 
effect is a consequence of the interference of partial waves with different angles of 
incidence. Creating a plane wave thus does not suffice if one wants to measure the 
GHS. But by superposing two plane waves generated according to (jSJ) with different 
4>(j) functions corresponding to a small difference in their propagation directions leads 
to a beam just like the one assumed in the derivation of the Artmann result: 

i/j(x, y) = if>i(x, y) + ip2(%, y) ~ exp (ikpix) + exp (ikp 2 x) , (6) 

with p m = sinx m ~ S mjX , where S mjX is the x component of the Poynting vector of 
partial wave m (with m = 1, 2). 
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As both partial waves of the beam §5§ are calculated using the same spherical wave 
components ipj, 

N/2 N/2 
i=-jV/2 J=-JV/2 

the sum is just given by 

N/2 

i/j = fa + fa = E V»i (e i0l(i) + e i</,2(j) ) • (8) 

j=-N/2 

In the experiment, measurements from N = 18 antennas with positions Xj = 
xa + jd/V2, Uj = Va — jd/V% with xa = 40 cm, i/a = 5 cm and d = 0.5 cm are 
superposed; the center of the coordinate system is the center of the Teflon plate (see 
figure [2]). Because of the high resolution needed, it is not practical to measure the 
electric field on the whole Teflon plate, as it would take several years for the 18-antenna 
array. Instead, they are only measured on parts of it as indicated in figure El One such 
part lies on the incoming beam and one on the outgoing beam; the averages in (CQ) are 
performed over these fields. The third measured part is directly on the boundary of the 
Teflon plate. While it is not necessary for the extraction of the GHS, the data measured 
here allow the direct observation of the reflection process. Figure [2] shows the measured 
real parts Re(^) of the beams generated according to ([5]) with <pj = at a frequency of 
v — 15 GHz (corresponding to a vacuum wavelength A ~ 2 cm). The boundary of the 
Teflon plate is shown, as well as the position of the j = antenna in the middle of the 
antenna array. It is clear from figure [2] that both the incoming and the outgoing beams 
have nearly plane wave fronts and travel at an angle of approximately 45 degrees to the 
plate boundary. 

4. Results 

4-1. Generated "plane wave" beams 

Details of the beams in the different measured parts of the Teflon plate are shown 
individually in figure [3] together with their Poynting vectors. The positions r; n and r out 
given by 

/33.5\ /-19.0\ 
r in = ^ L(J j cm, r out = y g g ) cm (9) 

for the Poynting vectors in this figure. 

The spatial resolution of 5 mm is clearly sufficient to see the structure of the 
incoming and outgoing beams, which in fact approximate plane waves. Their respective 
Poynting vectors show a well-defined propagation direction. The incoming angle 
calculated from the averaged Poynting vector is x ~ 47 degrees, the respective outgoing 
angle is x ~ 45 degrees, which further supports the claim that, in fact, a plane wave 
travelling at an angle of 45 degrees has been created. The Poynting vectors of the beam 
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Figure 2. Measured real parts Re(V>) of the beams at v = 15 GHz. The plate is 
indicated as the gray region; the position (xa,Ua) of the middle (j = 0) antenna in 
the antenna array is shown as well. 



at the boundary reveal - perhaps surprisingly given that so little can be seen in Keip 
itself - that there are incoming and outgoing parts of the wave at the boundary. A part 
of it is reflected at the boundary, but another part penetrates outside the Teflon plate. 
The penetration depth seems to be a bit larger than the wavelength of 2 cm. 

The generation of plane waves with one propagation direction thus works, at least 
for high v values {y > 10 GHz; below that value, the wave fronts are less well defined 
and the extraction of Poynting vectors is thus not possible with high accuracy) . 



4.2. Extracted GHS 

Two plane waves generated according to the scheme discussed in the previous section 
can now be superposed. The GHS of the resulting beam can be extracted and compared 
to the Artmann result. Table [1] shows the different phase functions (fii(j), foti) used 
in this section and the angle of incidence x °f the beam according to (jBJ). The choice 
of the 4>(j) functions is of course somewhat arbitrary; here, they are chosen such that 
the j dependence is simple, the difference between <fii(j) and 02 (j) is small (as this is 
the approximation in the Artmann result), and such that a sufficiently broad range of 
angles of incidence results. 
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Figure 3. Real parts Re(i/>) of the beams (color scale) and Poynting vectors (black 
arrows) in the three measured parts of the Teflon plate as shown in figure [2] (a) 
outgoing beam, (b) incoming beam, (c) beam at the plate boundary. The average 
Poynting vectors of the incoming and outgoing beams are shown as red arrows. 



Table 1. Phase functions </>i,2(j) and the resulting incoming angles \ m degrees as 
extracted from the experimental data. 



Mj) 8j/9 2j/3 j/2 j/10 -j/2 -8j/9 

Mf) J 3j/4 j/3 j/4 -j/4 -j/3 -j 

X (degrees) 22.5 27.9 35.1 41.4 44.9 47.6 55.9 68.9 



The extracted GHS for the incoming angles given in tabled] is shown in figure H] for 
v = 10 GHz (A ~ 3 cm) and v — 15 GHz (A ~ 2 cm) together with the Artmann result 
( flOl) . which for TM polarization is given by 



nk sinx a/sui 2 x - 1/n 2 ' 

where k = 2tt/X is the wave number. kAs is thus independent of the wavelength. The 
error bars are the errors in (S) in ( out ), which in turn, as (S) in ( out ) is an average value, are 
given by the standard deviation. The error 5 As in As is then given by error propagation; 
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Figure 4. GHS kAs as extracted from the measured data for v = 15 GHz (red dots) 
and v = 10 GHz (blue triangles). The black solid line is the Artmann result, the black 
dashed line marks the critical angle for total internal reflection for n = 1.44. 



the errors are all approximately 5As/As ~ 20 % (incoming angles above the critical 
angle) and 5As/As ~ 50 % (incoming angles below critical angle). As the GHS below 
the critical angle is small, larger relative errors are expected. 

We see clear signatures of the GHS and the main features of the Artmann result 
(zero GHS below the critical angle, maximum GHS at the critical angle, approximately 
constant GHS above the critical angle, independence of kAs of k) are well captured. At 
v = 10 GHz, the GHS is systematically higher than the Artmann result predicts. This 
is not surprising, as deviations due to the finite width of our generated beams can be 
expected to become larger for smaller frequencies. 

In fact, the agreement of the measured GHS with theoretical predictions gets even 
better if the finite width of the beams is taken into account. The individual beams are 
not completely plane waves, as they have a width fixed by the width of the antenna 
array. If one approximates them as Gaussian (which is justified as Lai et al. [2] have 
shown that the precise beam profile has little influence on the resulting GHS) with 
a beam width given by the experimentally extracted value a ~ 3 cm, the resulting 
GHS As can be calculated numerically as described in [13] . Figure [5] shows the results 
together with the experimental data for v = 10 GHz and v — 15 GHz. The deviations 
above the critical angle are indeed explained very well by a Gaussian (non-plane wave) 
beam profile. 

Overall, we find convincing agreement between the GHS extracted from our 
experimental data and the theoretical prediction of the Artmann result. The deviations 
of beams with smaller frequencies from the Artmann result are due to the finite width 
of our generated beams, as a numerical calculation of the GHS using a Gaussian beam 
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Figure 5. GHS kAs extracted from the measured data at v = 15 GHz (red dots) and 
v = 10 GHz (blue triangles). For comparison, the calculated GHS for a Gaussian wave 
packet with the corresponding frequency is shown as well (red solid and blue dotted 
curves). The black dashed line marks the critical angle for total internal reflection for 
ri = 1.44. 



profile with a beam width corresponding to the one found in our experiments shows. 
5. Conclusions and outlook 

We presented measurements of the Goos-Hanchen shift of microwave beams at a 
dielectric interface (Teflon to air). The beams were generated by superposing the 
spherical waves produced by 18 microwave antennas, which leads to beams with 
approximately plane wave fronts. The GHS is extracted from the Poynting vectors 
of the incoming and reflected beams, which can be found from measurements of the 
electric field of the beams. Our results agree well with the classical result of Artmann. 

Our experiment demonstrates that GHS measurements in dielectric microwave 
cavities are possible and meaningful. A similar setup using not a rectangular Teflon 
plate, but one with curved boundaries, could be used to study the effects of boundary 
curvature on the GHS. When the GHS is incorporated into ray models, it is usually 
assumed that the curvature effects are small, and the GHS as calculated at a planar 
interface can be used. Microwave measurements could give insight into the question if 
and when this is possible, and what deviations can be expected. Another interesting 
experiment would be to excite modes in a microwave cavity and, by measuring 
electric fields, directly observe effects of the GHS which have been predicted in optical 
micro cavities, e.g. the formation of new modes in an elliptical cavity [11] or the phase- 
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space shift of periodic orbits [T3] . 
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